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ABSTRACT 

The possibility that mass ejected during Asymptotic Giant Branch (AGB) stellar 
evolution phases falls back towards the star has been suggested in applications rang¬ 
ing from the formation of accretion disks to the powering of late-thermal pulses. In 
this paper, we seek to explicate the properties of fallback flow trajectories from mass- 
loss events. We focus on a transient phase of mass ejection with sub-escape speeds, 
followed by a phase of a typical AGB wind. We solve the problem using both hydro- 
dynamic simulations and a simplified one-dimensional analytic model that matches 
the simulations. For a given set of initial wind characteristics, we find a critical shell 
velocity that distinguishes between ’’shell fallback” and ’’shell escape”. We discuss the 
relevance of our results for both single and binary AGB stars. In particular, we discuss 
how our results help to frame further studies of fallback as a mechanism for forming 
the substantial population of observed post-AGB stars with dusty disks. 

Key words: stars: AGB and post-AGB — stars: winds, outflows — methods: nu¬ 
merical 


1 INTRODUCTION 

The possibility that some material ejected during the AGB 
mass loss stage falls back onto the star has been the sub¬ 
ject of a number of studies, motivated by both theoretical 


and observational considerations ( 

Bujarrabal Alcolea & Neri 

1998 

|Soker||20011 |Van Winckel et al.| 

1998 Van Winckel 

1999 

[Waters, Trams & Waelkens|1992 

Zijlstra et al.|2001 1. 


The discovery of a population of post-AGB stars host¬ 
ing dusty, long-lived disks is a puzzle for standard stellar 
evolution theory as it is not clear how these disks originate 
i Glayton et al.|2014 Su et al.|2007 Venn, K. A. et al.|2014 |. 
Given that some occur in binaries, it has been assumed that 
most if not all, are the result of some form of mass transfer 


(Bujarrabal et al. 2005 Dermine et al. 2012 

Hinkle et al. 

2007 Malek & Cami||2014 Van Winckel et al. 

20091. 


Theoretically, there is diverse interest in AGB mass-loss 
fallback. Asking when the return of matter can significantly 
influence post-AGB evolution, Soker| (20011 considered the 


forces acting on back-flowing material in order to set lim¬ 
its on when such fallback can occur (though no trajecto¬ 
ries were calculated). In order for back flow to play a role 


* E-mail: zchen25@ur.rochester.edu 
f E-mail: afrank@pas.rochester.edu 
J E-mail: blackman@pas.rochester.edu 
S E-mail: nordhaus@astro.rit.edu 


in post-AGB and Planetary Nebulae (PN) evolution, that 
study concluded that the fallback mass should be high, the 
ejected velocity should be low, and that the mass should be 
concentrated along the AGB equator. 

For a spherically symmetric case, thermal pulses or fuel 
limited relaxation oscillations may be a source of such mass 
loss ( [Van Horn et al.]|2003[ ) . In [Frankowski fc Soker| ( |2009[ ) , 
backflows were invoked as a means of creating very late ther¬ 


mal pulses (VLTP). Hajduk et al. (20071 also considered the 
possibility that mass accretion can induce a VLTP in the 
context of the old nova CK Vul where the mass is assumed 


to come from a companion. Frankowski & Soker (20091, how¬ 
ever, considered the accreted mass falling back from AGB 
wind material that had become part of a planetary nebula. 


In Soker (20081, the possibility for fallback to occur via 
a stellar ’’evanescent zone” was explored where, in addition 
to the escaping wind, gas parcels do not reach the escape 
velocity but rise slowly through the zone and then fall back. 
Wind and bound gas were found to exist simultaneously out 
to distances of ~ 100 AU. 

Backflow has also been proposed as a mechanism for 
slowing evolution along post-AGB evolutionary tracks ( |Zi-| 
jlstra et al.||2001| and for depleting post-AGB dust stars of 


refractory elements (e.g. Van Winckel et al. (19981). 


Explaining the observations of disks around post-AGB 
stars (discussed above) has, however, been one of the main 
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reasons for exploring studies of fallback during AGB mass 
loss. Such a suggestion was part of the motivation driv¬ 


ing the Soker (20011 paper. In addition, Akashi & Soker 
(20071 suggested that circumbinary rotating disks can also 


be formed from slow-AGB-wind material that is pushed back 
towards the star by wide jets emanating from one of the stars 


(i.e. collimated PPN or PN outflows Bujarrabal Alcolea & 

Neii|([l9Ml)). 


A natural means for generating fallback in the context 
of binary stars however may be common envelope (CE) evo¬ 
lution. Both simulations and analytic models indicate that 
some fraction of the material expelled from the envelope 
may not escape the binary system dNordhaus & Blackman 


2006| [^rdhaus et al.||2007[ [Kashi fc Soker||2011| [Ricker fc 


Tamm 

2012 

Passy et al. 

2012 


therefore, expected to fall back after the CE mass ejection 
event. This can occur even when the CE event involves a 
low-mass object such as a massive planet being swallowed 


by an AGB star( Nordhaus & Blackman 20061. Planetary 


companions to main-sequence stars are plentiful and those 
orbiting within ~ 10 AU can be expected to plunge into 


their host stars during the giant branches (Nordhaus et al. 


2010 2011 Nordhaus & Spiegel 20131. In such cases, ma¬ 


terial would be ejected over short timescales (a shell) with 
sub-escape velocities (14 < 14sc). Thus, it is possible that 
circumbinary disks are the result of CE ejections in which 
some fraction of the envelope is ejected at sub-escape speeds 
which may drop back to be shaped into a disk. 

Given the interest in fallback on the AGB, in this pa¬ 
per we attempt to gain some insight into the dynamics of 
fallback shells by considering only the radial trajectories of 
the expelled gas. Of particular interest is the response of 
a sub-escape velocity shell to internal driving by the AGB 
wind which would likely follow the shell ejection. Note that 
we do not explicitly specify or model the origin of the mass- 
loss history as our goal here is simply to study the dynamics 
of the shell after its initial ejection from the star. In Sect. 
2, we introduce our hydrodynamic model and methods to 
solve the problem. In Sect. 3.1, we present the results of full 
hydrodynamic simulations. In Sect. 3.2 we present a semi- 
analytic model and compare it with the simulation results as 
well as present a zeroth order estimate of the critical velocity 
of the shell that distinguish collapse or escape trajectories. 


2 METHOD AND MODEL 

We use the AstroBEAR adaptive mesh refinement (AMR) 
code to solve the hydrodynamic equations (Cunningham et 


al. 2009 Carroll-Nellenback et al. 20131. AMR is a com- 


r\dimension 


putational technique that divides a computational cell into 
child cells when the physical condition suggests 
to (usually high gradient). 


|+V.(,v) = 0 
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Figure 1. First panel shows the boundary condition of the star 
when it emits stellar wind. Second and third panels show time 
dependence of density and velocity. 


where n = p/rumean is the particle number density and 
rrimean = 1.3 ruH is the mean particle mass. 

We begin with an AGB star at the origin and allow it fill 
the grid with a pre-shell wind. The mass of the star is 1 Mq. 
The escape velocity at 1 AU from the center of the star is 
Vesc = 15.43 km/s. We account for radiation pressure within 
the gas and dust by using an the effective stellar gravity of 
aGM^/C with a — 0.134. However, we do not specify the 
dust species and evolution in the simulation, that is, we only 
assume a steady state distribution of dust. The whole region 
can be optically thin which will be examined. 


frad — 


p{r)KtotaiS{r) _ (1 - a)GMQp(r) 


( 4 ) 


therefore, 


rSOAU rSOAl 

/ p{r)Ktotaidr = / 

JlAU JlAU 


47r(l — a)GM p{r)c 


dr 

( 5 ) 

where p{r) is the density of wind at radius of r, Ktotai is the 
total mass weighted opacity. Sir) is the radiation intensity, 
L is the bolometric luminosity of the AGB star and r is the 


P = nkbT 
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optical depth. In our simulation, the pre-shell wind is driven 
into the grid at ro = 1 AU with velocity vo = 20.39 kmjs. 
The wind has an injection density of po = 7.5 x 10~^^g/cnT 
yielding a mass loss rate of Mo = 6.84 x 10“^ MQyr~^. By 
mass conservation law, 

poVo{AU)^ = p{r)v{r)F (6) 


of AMR. Thus the ’’effective” resolution of our simulations 
is (3840)^ 


3 RESULTS 

3.1 Simulations from AstroBEAR 


we can substitute equation (§ into equation |5f and for com¬ 
putational convenience, eliminate vo and v(r) on both side 
since they are on the same order of magnitude. For a typical 
L = 1500 Lq AGB star, we can get that r < 1. Obviously, if 
we lower all the wind density (as well as the mass loss rate) 
by a factor of several, the whole region would be optically 
thin and the hydrodynamic property would not change since 
equations all scale linearly with the density. 

The simulation is isothermal and the temperature is 
100 K. This yields an initial Mach number of Mach = 22.4. 
The boundary of the star is a sphere and the inner condition 
do not change upon the variation of outer condition. We 
divide the inner region of the star into two parts, a solid 
sphere and a spherically concentric shell (Shown in the first 
panel in Figurel.). The solid sphere has time dependent but 
uniform density and zero velocity while the concentric shell 
has the same time dependent uniform density but linearly 
increasing velocity from the inner side to the outer side. The 
inner side velocity is zero and the outer side velocity is the 
wind velocity. R and z axes are both reflective boundaries 
while top and right edge of the computational box only allow 
outflow, so any negative velocity at the boundary is set to 
be zero. The free fall time at the boundary is 80 yrs and the 
sound speed is 0.90 kmjs. Given that the initial ambient 
density is extremely low and that constant wind will fill the 
grid before the dense shell outflow, the boundary condition 
will have tiny effect on the simulation. 

After the wind fills the grid we drive a dense shell 
from the star. The shell episode lasts from ti = 2.32 yr 
to t 2 = 2.90 yr. During the episode, the density rises to 
1.5 X 10“^® gjavA. The velocity drops to either = 
14.27 kvajs or Vuigh = 14.68 km/s for the two cases we 
present here. As we will see, these two cases bracket a 
bifurcation in the evolution of the ejected shell. Both of 
these cases have Vs < Vesc and in both cases approximately 
M = 5.6 X 10“® Mq is ejected in the shell. For computa¬ 
tional convenience, we approximate the time dependence of 
density and velocity during shell ejections as step functions. 

After the shell ejection terminates, we resume the 
steady wind (which we subsequently refer to as a ’’post-shell 
wind”). The second and the third panels in Figurel show the 
density and velocity history for our simulations. The figures 
show that in both simulations, the velocity during the shell 
ejection is sub-escape at the stellar boundary where mass 
flows into the grid. 

Because AstroBEAR does not have 1-D capacities, we 
run our code in (r, z) cylindrical symmetry (2.5U) covering 
one quarter of the total plane and average over the spherical 
angle 9. Thus we impose reflecting boundaries at r = 0 and 
z = 0. The top and right boundaries use outflow boundary 
conditions allowing material to leave the grid. The winds 
and shell are injected into the grid via a spherical ’’surface” 
boundary condition at ro- The computational domain is 24 x 
24 AU^ and consists of a 240 x 240 base grids plus 4 levels 


In what follows we present two examples from a larger suite 
of simulations chosen to bracket the transition from escape 
to fail-back modes of the ejected shell. In Figure 2 and 3, we 
present ID profiles of density and velocity at four different 
times in the evolution of the simulations. We take the profiles 
along the 45-degree line in the r — z plane. We also plot the 
local escape velocity for reference. 

The four panels in Figure2 show frames from the 
simulation with the lower shell ejection velocity Vaheii = 
14.27 kmjs (note with time increases from top to bottom). 
The simulation begins with the steady wind filling the grid 
then the pulse of high density, sub-escape material (the 
shell) is launched at 2.32 yr. In panel (a) we see the density 
and velocity profiles at at t = 5.23 yr. Note that the steep 
rise in velocity between the origin and r = 1 AU represents 
the inflow boundary where the stellar flow is injected into 
the grid. In this region we also see that the density is held 
constant. 

In the first panel between 1 AU and 4 AU, the density 
of the steady wind ejected after the mass loss ’’pulse” drops 
continuously due to the geometric expansion of the gas. The 
velocity also falls due to deceleration by the gravity of the 
AGB star. Thus since p oc Mn,/Vw{r)r‘^, the wind density in 
this region falls off more slowly that p(r) oc r~^. 

At r = 4 AU, the steady post-pulse wind encounters 
the inner edge of the shell ejected during the brief pulse 
high mass loss. From this shock interface to approximately 
r = 6 AU we see shell material. Note that the velocity in 
this region is below the local escape velocity of the AGB 
star. Between r = 6 AU and r = 9 AU we see a steady drop 
in the density and steady increase in the velocity. This re¬ 
gion represents a rarefaction wave running from the denser 
but slower moving shell material and the supersonic, supra- 
escape pre-shell wind ejected before the pulse. The most 
important point to note in panel (a) of figure 2 is the pres¬ 
ence of the shock at r = 4 AU. This shock wave represents 
a transfer of momentum from the post-pulse wind into the 
shell which causes varying degrees of shell acceleration de¬ 
pending, on the relative momenta of the wind and shell. 

At t = 11.04 yr, we see the dense shell is has expanded 
outward. The shock bounding the shell and the post-shell 
wind is now at r = 6 AU. The boundary between the shell 
and pre-shell wind has expanded off the grid. Note however 
that the velocity of the shell is just above u = 0. Thus the 
shell has been almost entirely decelerated and is reaching its 
maximum radial excursion. 

By t = 16.86 yr, the velocity of the shell has fallen be¬ 
low II = 0. At this point the shell has begun to fall back 
towards the star. By t = 22.09 yr, the shock boundary be¬ 
tween the wind and the shell is at r = 3 AU as the shell 
continues to collapse back towards the star. Note that for 
all positions above this radius we see n < 0 which implies all 
material at these larger radii is moving radially inward as 
well. This makes sense because of the rarefaction between 
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Figure 2. Time evolution of shell with initial shell velocity Ujou, = 
14.27 km/s. The shell finally falls back. 


the shell material and the pre-pulse wind. Since the right 
and top boundary is made to only allow outflow, it must 
be different from a self-consistent boundary. However, our 
simulation zone is large enough to keep the region we are 
interested unaffected. If the gas flows in from the boundary 
right after the velocity become zero, an error - induced by 
the artificial boundary condition - will travel at the sound 
speed which is 0.90 fem/s. However, the error does not have 
enough time to travel from the boundary (24 AU) to 12 AU 
in just 22 yrs. So the fluid motion in 12 x 12 AU^ is not 
adversely affected by the boundary. 

Finally, note that the sound speed becomes greater than 
shell speed as it decelerates (cs > v). During this period 
we have a subsonic flow {Mach < 1) where pressure force 
become important and we see an expansion of the shell width 
with very slow bulk radial motion. 

In Figure 3, we show radial profiles of density and ve¬ 
locity for the simulation with a higher initial shell velocity. 
In the first panel, taken at t = 10.46 yr, the shock boundary 
between post-shell wind and the shell is located at about 
r = 7.5 AU. Note the velocity of the shell is below the 
local escape velocity. As the post-shell wind continues to 
impart momentum to the shell, we observe a different evo¬ 
lution. In subsequent frames, the shell does not stall but 
is slowly driven outward. By t = 31.39 yr, the shell has 
reached r = 16 AU and its velocity is now just below the 
local escape value. By t = 41.85 yr the shell has been al¬ 
most entirely pushed off the grid and has achieved a speed 
in excess of the local escape velocity. Thus we find a shell 
ejected with Vshm = 14.27 km/s, will fall back onto the star. 
The shell with Vsheii = 14.68 km/s, which is also below the 
escape velocity, at the launch radius, will escapes. 


3.2 A One Dimensional Spherical Model 

In order to understand the transition between models with 
shells that escape under the action of the post-shell wind 
and those which fall back onto the star, we now derive a 
simplified analytical treatment of the problem. 

To compare our 2.5D simulation with a ID analytic 
model, we use fluid tracers in the simulations to track 
shell material. We use tracer mass weighted averages of the 
tracer’s radial extent as a measure of shell’s position vs time 

rs{t). 

fiOAU 3^ 

_ JlAf/ Ptracerl ai 

rSOAU 9 T 

j Ptracer'k' dv 

We assume that the initial velocity and density of the shell 
are Vs and ps respectively. For the shell initial mass we have: 

ms = A-Kr^ psVsAt, (8) 

where At is the time interval of injection. The post-shell 
wind has initial velocity vo and initial density po at the 
launch point. For a steady spherical wind, mass conservation 
implies: 

47rr^p„u™ = 47rropoUo (9) 

Thus the wind density will fall off as: p™ = (roPoUo)/(j'^Um). 

We calculate the evolution of the shell velocity Vs from 
the momentum equation it . We assume that pressure forces 
are negligible with || — VP|| <C || ^ II • This assumption 
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Figure 4. Results of trajectories of the shell from AstroBEAR 
and analytic model. Pulsation velocities are viow = 14.27 km/s 
and Vhigh = 14.68 km/s respectively. 


is valid at all times except during the brief stalling of the 
shell that initiates fallback. The momentum equation thus 
becomes 


dp^ 


+ V • (pvv) = 


GMp. 

——r. 


( 10 ) 


We must track the velocity of the wind gas before it impacts 
the dense shell as it is subject to deceleration due to gravity, 


Vw{r) 



2GM 

ro 


2GM 

r 


( 11 ) 



For the momentum of the shell we write a discrete form 
of the momentum equation as 

/ , \ \ \ \ , GMms , , 

ms{t + dt)Vs(t + at) = ms(t)'Vs{t) + dmw'Vm - ^—rdt 

( 12 ) 

where the second term on the right represents momentum 
added to the shell by the wind in a time dt and the third term 
represents deceleration due to gravity. More specifically, for 
the wind mass added to the shell we have, 

dm-u, = 4:TvF p{v.u,{t) — Vs{t))dt (13) 

so that 

ms{t + dt) = ms{t) + dniw (14) 


Figure 3. Time evolution of shell with initial shell velocity 
'>’high = 14.68 kvajs. The shell exceeds the escape velocity 


Using the evolution of the shell velocity: 


Vs(t + dt) 


ms{t)vs + rdt 

ms{t + dt) 


(15) 


Using Eqs. (9) and (10), we divide equation (8) by 
ms{t + dt) and solve it numerically to find the trajectory 
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of the shell r{t) for any set of initial conditions. In particu¬ 
lar we look for the critical value of the initial shell velocity Vc 
that distinguishes between fall back and escape. Evaluation 
of equation (15 I shows that for the inputs used in the simu¬ 


lations, the critical initial shell velocity is Vc = 14.50 kmjs. 
Figure 4 shows a comparison of the rs{t) for the analytic 
model and for simulations (i.e. vo < Vc and vo > Vc). The 
analytic model does quite a good job of recovering the be¬ 
havior seen in the simulations in both cases. 

We can also estimate Vc even more simply by first ap¬ 
proximating the turning point of the shell in the fall back 
case. Here assume that the mass of the shell is constant and 
neglect the effect of post-shell wind. This gives us a ballistic 
approximation for the maximum radius of the shell r^. We 
then assume that the ram pressure of the steady wind must 
be large enough to support the shell against gravity at this 
radius if the wind is to stop the shell from falling back. This 
yields 

GtTIs^^ dpram 


dt 


(16) 


where Pram is the the magnitude of the ram momentum of 
the wind given by 

Prarn = 47rpr^u^ = Atv porlvov,^ = Mvm (17) 


Note that Vw is a function of r from Eq. (111. Note also 


that rm is determined by vq according to our assumption of 
ballistic motion, namely it corresponds to the radius where 
the velocity vanishes. Then force balance gives 

-2roGM 


vgro - 2GM ’ 


(18) 


2GM 2GM 

- + -, 

ro rm 


Then using Eqs. ( [T8| ), ( [T^ and|^ in the form 
rus = inrlpoVcAt 


(19) 


( 20 ) 


in Eq. (161, we find Vc- The result is Vc = 14.77 kmjs which 


is only 1.9% higher than the numerical result 14.50 km/s 


4 SUMMARY AND DISCUSSION 

We have presented a simple model for the evolution of a 
shell ejected with sub-escape velocity during a brief pulse of 
enhanced mass loss in an AGB star. The fate of the ejected 
shell depends on the acceleration to escape velocity via the 
action of the continuing steady AGB wind that follows the 
launching of the shell. Our 2.5D isothermal hydrodynamic 
simulations are compared with a spherically-symmetric an¬ 
alytic model for the time evolution of the shell. Both the 
analytic model and simulations closely support in their pre¬ 
dictions for the value of the critical shell velocity Vc that de¬ 
lineates the bifurcation between escape and fail-back modes 
of shell evolution. Though in the presence of spherically sym¬ 
metric radial temperature gradient, the diffusion rate will be 
different but the qualitative result that we hnd-a bifurcation 
of fall back vs. escape modes would still be expected. 

Our results lay out the basic dynamics of fail-back shells 
and are relevant to any mechanism which ejects material 
with less than the local escape velocity. Our study also es¬ 
tablishes the first step along the path to study post-AGB 


disk formation via fail-back shells. Since AGB stars are slow 
rotators, we would expect that an ejected shell would not 
have enough angular momentum to establish a Keplerian 
disk at large radii. Thus ejected gas will either fall back onto 
the star or escape for isolated stars. In the case of a binary 
however, the material falling back onto the AGB star might 
gain enough angular momentum through gravitational in¬ 
teraction to form a Keplerian disk at large radii. This is a 
topic for future study. 
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